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Raman amplification in the coherent wavebreaking regime 

J. P. Farmer and A. Pukhov 

Heinrich Heine Universitdt, 40225 Diisseldorf, Germany 

In regimes far beyond the wavebreaking threshold of Raman amplification, we show that significant 
amplifcation can occur after the onset of wavebreaking, before phase mixing destroys the coherent 
coupling between pump, probe and plasma wave. Amplification in this regime is therefore a tran¬ 
sient effect, with the higher-efficiency “coherent wavebreaking” (CWB) regime accessed by using a 
short, intense probe. Parameter scans illustrate the marked difference in behaviour between below 
wavebreaking, in which the energy-transfer efficiency is high but total energy transfer is low, wave¬ 
breaking, in which efficiency is low, and CWB, in which moderate efficiencies allow the highest total 
energy transfer. 


I. INTRODUCTION 

Raman amplification in plasma has been suggested as 
a mechanism to allow the creation of ultrashort, ultrain¬ 
tense laser pulses, [1] which have applications across sci¬ 
ence and technology. Using plasma as a gain medium of¬ 
fers the advantage that, unlike solid-state media, it does 
not have a damage threshold, potentially reducing or re¬ 
moving the need for stretching and compression of the 
laser pulse, allowing either a reduction in size or an in¬ 
crease in power of next-generation laser systems. 

In the Raman interaction, a probe pulse interacts 
with a counterpropagating pump, driving a plasma wave 
through the ponderomotive force of their beat. If the 
pump and probe frequencies, uja and are chosen such 
that the probe is downshifted from the pump by the 
plasma frequency, the plasma wave is resonantly ex¬ 
cited and may grow to large amplitude. The resulting 
density perturbation acts as a moving Bragg grating, 
which acts to scatter and Doppler-shift the pump pulse 
into the probe, amplifying the latter. 

Although the plasma itself does not have a damage 
threshold, there remain several effects which act to limit 
the interaction, such as filamentation and parasitic spon¬ 
taneous backscatter of the pump [2]. Further, the high 
wavenumber and low phase velocity of the laser beat, 
which allows a large-amplitude plasma wave to be ex¬ 
cited, means that the plasma wave can break for even 
moderate laser intensities, <C 1x10^^ W cm“^ [3], lead¬ 
ing to strong damping [4], which will limit the coupling 
between pump and probe [5]. 

These effects likely contribute to the low efficiencies 
reported by experimental campaigns [6-9], the best to 
date being by Ren et a/., achieving 6.4% in a double-pass 
setup [10]. Previous theoretical and computational works 
have focussed on the optimal pump amplitude for efficient 
amplification [2, 5, 11, 12]. Although it is known that the 
highest efficiencies are achieved below the wavebreaking 
threshold, the necessarily low intensities would require 
large interaction volumes to scale to high power. In this 
work we therefore focus on the energy-transfer efficiency 
far beyond the wavebreaking threshold. 

It is interesting to note, however, that simulations car¬ 
ried out in this regime have yielded significantly differ¬ 


ent efficiencies, even for similar parameters, ranging from 
35% in the work of Trines et a/. [2], to less than 10% in 
that of Toroker et a/. [11]. These apparently contradic¬ 
tory results have been the subject of significant discus¬ 
sion [11, 12], although the main cause has not previously 
been identified. Clearly, reconciling the differences be¬ 
tween simulations is important if reliable comparison to 
experiment is to be made. 

In this work we identify a new process in the Raman 
interaction, coherent wavebreaking. In section II, we il¬ 
lustrate a strong dependence of the energy-transfer ef¬ 
ficiency on the probe duration, which explains the dif¬ 
ference between efficiencies observed in other works. A 
simple analytical model is developed in section III, al¬ 
lowing the relevant physical processes to the identified. 
Section IV discusses the applicability of the results to 
experiments, and conclusions are drawn in section V. 


II. INFLUENCE OF PROBE DURATION FAR 
BEYOND THE WAVEBREAKING LIMIT 

We make use of the Leap code [13], which is based 
on a laser-envelope particle-in-cell model. Since the Ra¬ 
man interaction is predominantly planar due to the short 
wavelength of the excited plasma wave and the relatively 
large interaction cross section, we here limit ourselves 
to one-dimensional (ID) simulations. In addition to a 
low computational overhead, this geometry has the fur¬ 
ther advantage that the plasma response modelled by 
the Leap code is exact [14]. Although multidimensional 
simulations are certainly important for direct compari¬ 
son to experiments, the fundamental processes are often 
obfuscated, as different regimes of amplification may be 
in effect at different radii, discussed in more detail in 
section IV. We therefore make use of ID simulations in 
order to characterise the coherent-wavebreaking regime 
of interest here. 

Figure 1 shows the evolution of the probe pulse for dif¬ 
ferent initial probe durations. Gaussian pulses are used, 
with FWHM-intensity durations of 17, 50 and 83 fs. A 
flat-top, 800 nm pump with duration 30 ps and inten¬ 
sity 1x10^^ W cm“^ is used to amplify a probe with 
initial intensity 1x10^^ W cm“^ in a plasma of density 


2 


<u 

"O 


Q. 

E 

< 



Position (fs) 


E 

u 



-M 

'to 

c 

<u 


6x10^^ 

3x10^^ 


0 



50 


25 


0 


m 

(D 


Figure 1. (Color online) Plots showing evolution of the probe pulse for three different initial durations, a-d) Snapshots of the 
amplitude at 0.5, 5.5, 10.5 and 15.5 ps. e) Peak probe intensity against time, f) Total probe energy against time. 


4.4 X10^^ cm“^ with an initial temperature of 10 eV. This 
corresponds to uja /= 20 , uJhl^p = 19, and a pump in¬ 
tensity 30 times above the threshold for wavebreaking 
to occur [3]. These are the same parameters as used by 
Toroker et al. [11] in their investigation of the strong¬ 
wavebreaking regime, with the longest pulse used here 
equivalent to the probe used in that work. Although the 
longer pump length used here (30 ps compared to r\j 1 ps) 
would result in higher plasma temperatures, we retain the 
10 eV temperature to allow direct comparison to those 
results. This longer pump allows comparison with the 
results of Trines et al 

Figures la-d show snaptshots of the probe profiles at 
different times during the interaction, given as the re¬ 
duced vector potential eEjmcuJi), with E the electric field 
amplitude, c the vacuum speed of light, and —e and m 
the electron charge and mass. It can be seen that in all 
cases the probe is amplified and compressed. However, 
as highlighted in Fig. le, the greatest amplification is 
achieved for the shortest pulse, which reaches an intensity 
of 6 xlO^^ W cm The reduction in peak intensity at 
the end of the simulation is due to dispersion as the probe 
continues to propagate through plasma after the end of 
the pump at 15.5 ps. Energies were calculated by scaling 
to three dimensions, taking a pulse radius (e-squared- 
folding distance) of 595 /rm, equal to that used by Trines 
et al. [ 2 ], and are shown in Fig. If. While the 17 fs probe 
is initially the least energetic and is compressed to 9.7 fs 
during amplification, it goes on to exceed the energy of 
the longer pulses. The total pump energy for these pa¬ 
rameters is 166 J, giving pump-to-probe energy-transfer 
efficiencies over the full interaction for the three probe 
durations of, from shortest to longest, 28, 21 and 14%. 

We note from Fig. If that the growth of the probe en¬ 
ergy is better than linear, with the efficiency increasing 
over the interaction as the probe becomes more intense. 
The gain of 23 J over a 30 ps pump observed for the 83 fs 


probe is therefore consistent with the results of Toroker et 
a/., which for this pulse width correspond to an energy- 
transfer of ^0.5 J over ^1 ps pump. Trines et al. used 
similar parameters to the 50 fs probe used here, but with 
cold plasma and an on-axis intensity for pump and probe 
double that used here. However, even with those factors 
taken into account, the amplification here remains some¬ 
what lower than observed in that work. This disparity 
could be a result of the initial pulse shape - using a pulse 
profile with a finite support, for example truncating a 
Gaussian at 2 cr, increases the amplification, giving com¬ 
parable results. 

The improved amplification observed for shorter probe 
pulses and the sensitivity to the pulse shape are linked 
to self-steepening, which causes the peak of the amplified 
pulse to move forward, away from the initial maximum. 
The effective seed amplitude is therefore reduced, as a 
point on the leading edge of the initial pulse becomes the 
seed for amplification, rather than than the peak. 

Self-steepening in the Raman process can arise due 
to many different effects, such as ponderomotive non¬ 
linearity [1], pump depletion [3] and the use of a chirped 
pump [15], while the presence of a pre-pulse can lead to 
superluminous precursors [16]. The self-steepening ob¬ 
served here is a result of wave breaking, as can be seen 
in Fig. 2 . Snapshots of the interaction 1.2 ps into the 
simulation, before significant amplification occurs, show 
the pump, probe, and plasma-wave amplitude, the elec¬ 
tron phase space and the electron density. The plasma 
wave amplitude is here given as the absolute value of 
the coupling susceptibility [13] normalised to the square 
of the plasma frequency, \ijj\/uj‘^ = | ((n/qno) |, 

where (pa^b are the vacuum phases of the pump and probe, 
n and no are the local and equilibrium plasma electron 
densities, and 7 the plasma-electron Lorentz factor. This 
value represents the coupling between pump and probe. 
From linear theory [4], wavebreaking occurs when the 
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Figure 2. (Color online) Snapshot of the Raman interaction at time 1.2 ps for initial probe durations of (from left to right) 
17, 50 and 83 fs. a) The pump and probe amplitudes and the effective coupling between the two. b-d) plots of the electron 
phase space on different scales, showing b) the long-wavelength oscillation excited by the breaking plasma wave, c) the onset 
of wavebreaking itself (ubr ~ 0.025c), and d) the electron phase space at the probe peak, e) shows the corresponding electron 
density at the probe peak. 


electron velocity exceeds the phase velocity of the wave, 
Vhr/c ~ ujpl2uja = 0.025, corresponding to \ip\/ujp = 0.5. 

For the longer pulses, wavebreaking occurs on the lead¬ 
ing edge of the probe. As a result, phase mixing of the 
plasma wave has already set in by the time the probe 
peak is reached, as seen in Fig. 2d. The resulting density 
perturbation, shown in Fig. 2e, loses its periodic struc¬ 
ture and no longer efficiently scatters the pump into the 
probe, lowering amplification. 

Although the electrons periodically rephase as they 
make synchrotron-like oscillations in the broken wave, 
the shift in resonance means that these plasma echoes 
have the wrong phase to coherently scatter the pump 
into the probe. Amplification therefore depends on the 
coupling susceptibility prior to phase mixing. For the 
shortest pulse in Fig. 2, the first peak in the coupling sus¬ 
ceptibility almost coincides with the peak of the probe, 
maximising the energy transfer. 

If the probe is too short, however, peak amplification 
may occur behind the probe peak, again lowering energy 
transfer. This is in fact the case for the 17 fs probe for 
time < 1.2 ps, giving a lower energy-transfer rate than 
observed for the 83 fs probe. However, as the probe is 
amplified, the point at which wavebreaking occurs moves 


forwards. This causes the point of peak coupling to 
advance, leading to improved overlap with the probe, 
increasing energy transfer. From 1.2 ps onwards, the 
energy-transfer rate is highest for the 17 fs probe, which 
goes on to become the most energetic after ^ 5.5 ps. 

In addition to the short-wavelength excitation, 
^ (TTc/cJa), driven by the beat of the pump and probe, 
it is interesting to note that there is a long-wavelength 
excitation, ^ (27rc/cc;p), as seen in Fig. 2b. This is not 
a wake driven by the probe, and is not observed below 
the wavebreaking threshold, and is instead driven by the 
breaking of the short-wavelength wave. From ID cold- 
plasma theory, below the wave-breaking threshold, indi¬ 
vidual electron charge sheets do not cross. The force act¬ 
ing on a sheet therefore depends only on the smeared-out 
ionic background, leading to simple harmonic motion[4]. 
When the wave breaks, charge sheets cross, and so ex¬ 
perience an additional force due to the electron-charge 
imbalance. Since the point at which the wave breaks 
travels with the probe, the associated plasma excitation 
has a phase velocity ~ c, resulting in a long-wavelength 
excitation. 

A similar effect has been identified in low-density 
plasma, in which the excitation arises from ponderomo- 
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tive nonlinearity, and has been suggested as a method 
to allow controlled electron acceleration with relatively 
low-intensity pulses. [17] The effect observed here has the 
potential advantage that higher density plasma may be 
used, allowing a larger accelerating field. Further, there 
is no constraint on the probe duration, making the effect 
more widely accessible; for a sufficiently intense pump, 
the only requirement is that the integrated probe ampli¬ 
tude is sufficiently high for wavebreaking to occur. 


III. ANALYTICAL MODEL 


The decrease in the effective probe amplitude caused 
by self-steepening is comparable to the effect of shadow¬ 
ing [18], in which the leading edge of the probe is prefer¬ 
entially amplified, lowering the effective seeding power. 
The effect here, however, is distinct, as wavebreaking 
changes the interaction, preventing the 7r-pulse solution 
which arises in the pump-depletion regime [3]. To il¬ 
lustrate this, we consider the coupled three-wave equa¬ 
tions widely used in the study of parametric amplifica¬ 
tion. Phenomenological treatments for wavebreaking in 
the three-wave model have been investigated in other 
works [6, 19], as have three-wave models incorporating 
damping and frequency shifts calculated from a nonlin¬ 
ear density distribution function [20]. Rather than at¬ 
tempting to derive the exact behaviour of the system in 
regimes far above the wavebreaking threshold, we find 
that using only simple assumptions we can recover the 
qualitative behaviour observed in Fig. 1. We note that, 
although dispersion plays a role towards the end of the 
simulation in Fig. 1, amplification remains the dominant 
process until the end of the pump-probe interaction. We 
therefore limit ourselves to the linearised, dispersion-free 
equations for pump, probe and plasma wave, to obtain: 


{dt -dz)a= 



(d^ + n);p = !!^a*b. ( 1 ) 

Here a and b are the envelopes of the pump and probe, 
respectively, which satisfy a = 91e ((ae*"^® iT), 

with a the reduced vector potential, 0 a ,6 = ^a,b{t 
the carrier phases of the pump and probe, and u the 
polarisation vector, u = {x ^ iy)jyj2 for circularly 
polarised light. We retain the coupling susceptibility, 
0 = !SQ^m) which is related to the nor¬ 

malised electric field used in, e.g., [3], by 0 = 2iuJa^pf'^- 
The functional f2[x,t,'0] allows the effects of wavebreak¬ 
ing to be taken into account. Assuming an initially 
cold plasma, valid for the low temperatures considered 
here, we expect no infiuence from wavebreaking below 


the threshold ^br = O.Scj^, which yields: 


f2[x, t, 0] 


0 , 




max 

max 


x,t') ) < ^br 





( 2 ) 


where u and S are the damping rate and shift in the res¬ 
onant frequency of the plasma wave due to wavebreaking 
(in a full treatment, these quantities will certainly vary 
in time). 

Assuming the probe is sufficiently intense for wave¬ 
breaking to occur (readily satisfied for high pump in¬ 
tensities), the analytical 7r-pulse solution for the pump- 
depletion regime will remain valid only for the leading 
edge of the interaction, before the wave breaks. Sub¬ 
stituting the wavebreaking threshold for 0 into these 
equations, we obtain the pump amplitude at wavebreak¬ 
ing, ttbr = (^0 — (co’p/2Co’a)^)^/^. Wc UOtC that ttbr 

is independent of the initial probe amplitude. More¬ 
over, we find that the pump depletion at wavebreaking, 
^ aQ — = (co’p/2co’a)^, is independent of both pump 

and probe amplitudes, and corresponds to the threshold 
pump intensity for wavebreaking to occur. 

The energy density of the plasma wave at wavebreak¬ 
ing depends only on its frequency and wavenumber. Since 
the pump depletion is fixed, it follows that the intensity 
increase of the probe up to wavebreaking is also indepen¬ 
dent of the laser amplitudes. Either the probe growth 
rate is high, in which case wavebreaking occurs rapidly, 
or the growth rate is low, and wavebreaking occurs pro¬ 
portionally later. The energy gain can only be modified 
by changing the laser and plasma frequencies, which al¬ 
ters the energy partition and the wavebreaking threshold, 
or by changing the interaction volume, i.e. a longer pump 
or wider-diameter beams. 

Therefore, in order to improve amplification, we must 
consider the interaction after wavebreaking. As the 
plasma electrons move with finite velocity, 0 must be 
continuous in time, and so Q must remain finite. Since 
the damping and the shift in resonance arise due to phase 
mixing within the broken wave, we expect the character¬ 
istic decay time of the wave to be of the order l/ujp. In 
the finite time for coupling to break down after wave¬ 
breaking, we therefore find the total energy transfer will 
be greater for larger pump and probe amplitudes. 

Put simply, the energy transfer prior to wavebreaking 
is independent of the laser amplitudes, and so efficient 
amplification relies on maximising the energy transfer af¬ 
ter the wave has broken. This is achieved by maximis¬ 
ing the probe amplitude in the period immediately after 
wavebreaking, before phase mixing destroys the coupling 
between pump and probe, which requires high-contrast 
probe pulses of duration ^ l/ujp. 

The “soft-wavebreaking” model discussed by Balakin 
[ 6 ], in which wavebreaking limits but does not reduce 
the probe growth, can be considered as equivalent to this 
model in the limit of very low plasma density, where the 
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Figure 3. (Color online) Pump-to-probe energy-transfer effi¬ 
ciency for varying pump intensity. The vertical dashed line 
shows the wavebreaking threshold. 

probe duration is ^ l/cjp. The model discussed here, 
however, remains valid for the higher plasma densities 
necessary for significant amplification. 

We therefore make a distinction between the wave¬ 
breaking regime [11], in which amplification of a long 
probe is significantly reduced by phase mixing after 
the plasma wave breaks, and the coherent wavebreaking 
(CWB) regime, observed for the shortest pulse in Fig. 1, 
in which the probe is sufficiently short that the probe 
peak is amplified, even after the wave breaks. The CWB 
regime is therefore characterised by much higher efficien¬ 
cies than the wavebreaking regime, and so the use of a 
high-quality probe pulse, with short duration and good 
contrast ratio, is vital for experimental campaigns. 

IV. APPLICABILITY TO EXPERIMENTS 

To illustrate this importance. Fig. 3 shows the pump- 
to-probe energy-transfer efficiency while varying the 
pump intensity for different initial probe durations. Only 
energy transferred to the probe within 500 fs the initial 
probe peak is considered. Other parameters are as for 
Fig. 1. Using a shorter pulse gives rise to higher efficien¬ 
cies over a wide range of pump intensities beyond the 
wavebreaking threshold. Despite this, better efficiencies 
are achieved by using a pump below the wavebreaking 
threshold. In this regime, longer probe pulses exhibit 
higher efficiencies, as the more energetic probe is better 
able to deplete the pump. 

Although below-wavebreaking amplification yields 
higher efficiencies, this regime is not necessarily prefer¬ 
able for experiments, as the large interaction volumes 
required for the same total pump energy may be techni¬ 
cally difficult to achieve. Physically, changing the pump 
intensity while keeping the probe constant, as is the case 
for Fig. 3, can be understood as the effect of stretch¬ 
ing the pump in time - the same total energy for a 
1x10^^ W cm“^ pump would therefore require a 45 cm 


interaction length. For limited interaction volumes, then, 
a higher intensity pump may be preferable in order to 
maximise energy transfer, despite the decrease in effi¬ 
ciency. With these considerations in mind, the use of a 
shorter probe pulse becomes important in order to access 
the the CWB regime. This regime has the additional ad¬ 
vantage that the duration of the amplified probe can also 
be significantly shorter than that achieved at lower pump 
intensities. 

Although the ^95% efficiency observed at low pump 
intensity is more than one order of magnitude better than 
the best experimental results to date, we note that no 
experiments have been carried out with similar param¬ 
eters - Ren et al. [10] used an initial probe intensity 
^10^^ Wcm“^, compared to the 1x10^^ Wcm“^ used 
here. We note, however, that the low efficiencies observed 
in recent campaigns at petawatt-scale facilities (Vulcan, 
PHELIX) may be due to the lack of a suitable probe - 
the shortest available pulses at such facilities are often 
^ 100 fs duration. 

Multidimensional simulations show that for beams 
with Gaussian transverse profiles, the optimal probe du¬ 
ration for a given pump intensity tends to be longer than 
that shown here, although the general trends are the 
same. This is due to the fact that at larger radii, the 
pump and probe intensities are lower, and so while the 
centre of the interaction may be far beyond wavebreak¬ 
ing, some parts of the interaction will be in the near¬ 
wavebreaking or below-wavebreaking regimes, for which 
longer probe durations are preferable. As each regime 
has different requirements for optimal efficiency, purpose- 
built Raman amplifiers would likely benefit from the use 
of fiat-top profile beams, as the optimal probe parameters 
are the same for all radii. 


V. CONCLUSIONS 

To conclude, we identify the coherent wavebreaking 
(CWB) regime of Raman amplification, in which signifi¬ 
cant amplification occurs after the onset of wavebreaking. 
The regime is accessed by using a short, intense probe 
pulse, which results in wavebreaking occurring close to 
the peak of the amplified probe, maximising the energy 
transfer from pump to probe in the time before signif¬ 
icant phase mixing occurs. This regime is of great im¬ 
portance to achieving high amplification, as experimental 
constraints act to limit the possible interaction volume, 
and as such limit the possible energy-transfer that can 
be attained below the wavebreaking threshold. 
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